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Diminution of the thermal conductivity is a crucial aspect in thermoelectric research. We report 
a systematic and significant reduction of the cross-plane thermal conductivity in a model system 
consisting of DC sputtered TiNiSn and HfNiSn half-Heusler superlattices. The reduction of n is 
measured by the Scj method and originates from phonon scattering at the internal interfaces. Heat 
transport in the superlattices is calculated based on Boltzmann transport theory, including a diffusive 
mismatch model for the phonons at the internal interfaces. Down to superlattice periodicity of 3 
nm the phonon spectrum mismatch between the superlattice components quantitatively explains 
the reduction of n. For very thin individual layers the interface model breaks down and the artificial 
crystal shows an enhanced n. We also present an enhanced ZT value for all investigated superlattices 
compared to the single TiNiSn and HfNiSn films. 

The ability to convert a temperature difference to elec¬ 
tricity as well as the possibility of both heating and cool¬ 
ing are very valuable properties of thermoelectric materi¬ 
als. Despite these promising features they are not widely 
used in industry as they posses low efficiency, described 
by a dimensionless figure of merit ZT = ^^T {S - See- 
beck coefficient, a - electrical conductivity, Ktot - thermal 
conductivity, T - absolute temperature). 

MNiSn half-Heusler (HH) materials, where M = (Ti, 

Hf, Zr) are considered to be materials having high poten¬ 
tial for thermoelectric (TE) applications due to Seebeck 
coefficients in the range of —200 /iV/K for bulk materi¬ 
als already at room temperature [1]. However, simulta¬ 
neously they exhibit thermal conductivities even as high 
as 10 W/(mK) preventing the achievement of a satisfac¬ 
tory value of XT, which is a measure of the material’s 
applicability [1]. Nowadays the scientific community is 
exploring all means to reduce the thermal conductivity of 
HH materials in order to enhance their ZT. Promising ap¬ 
proaches include: the introduction of grain boundaries by 
melt spinning or ball milling processes and further spark 
plasma sintering [2] , phase separation during the solidifi¬ 
cation of bulk materials [3-6], or the thin film and super¬ 
lattice (SL) approach [7, 8]. The latter was investigated 
by Venkatasubramanian et al. for the Bi 2 Te 3 /Sb 2 Te 3 SL 
system and resulted in the highest ever reported ZT value 
of 2.4 at room temperature [9, 10]. In contrast to tel- 
lurides, HH materials posses the peak of efficiency at high 
temperature (above 700 K) [11]. Moreover, the cost of 
bulk HHs is about one order of magnitude lower com¬ 
pared to the former compounds [12]. 
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The superlattice approach is a promising way to re¬ 
duce the thermal conductivity, however, the total ther¬ 
mal conductivity consists of an electronic and a phononic 
contribution (/^tot = + ^ph)- The first is directly re¬ 

lated to the electrical conductivity via the Wiedemann- 
Franz law (^ = const) [13]. Therefore, reduction of 
leads to a reduction of a and to a decrease of ZT. 
To ensure that and other electronic properties remain 
unchanged, isoelectronic elements are used as substitutes 
for the M element in the compounds so that mainly the 
phononic thermal conductivity is affected. In our model 
system SLs have constant total thicknesses, but variable 
SL period and therefore number of interfaces. However, 
SLs with ultrashort period consist only of a few atomic 
planes of the different materials. Thus they should bet¬ 
ter be considered not as layers of material 1 and material 
2 separated by interfaces, but rather as an artificial tai¬ 
lor made new material with a large crystallographic unit 
cell. 

TiNiSn/HfNiSn SLs having the same amount of both 
materials in one period were grown by DG magnetron 
sputtering. Details about sample preparation are pre¬ 
sented in the supplementary data. The cross-plane ther¬ 
mal conductivity of HH SLs was measured with the 3uj 
method [14-16]. The sample geometry and the measure¬ 
ment procedure was described in detail in our previous 
work [7] and in the supplementary materials. 

A crucial aspect of the approach pursued here is a 
strong sensitivity to the quality of the interfaces. X-ray 
diffraction (XRD), atomic force microscopy (AFM) and 
scanning transmission electron microscopy (STEM) were 
used to estimate interface quality. Assuming that the 
interface roughness transfers to the surface roughness, 
a root mean square (rms) roughness of 1.1 nm for a 
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FIG. 1. (a) STEM dark field image of a SL cross-section with a period of 4 nm, composed of TiNiSn (dark contrast) and 
HfNiSn (bright contrast) and surface topography of an approximately 800 nm thick SL with rms roughness of 1.1 nm (AFM 
image) in inset, (b) SAED of a circular region from the image (a), (c) high-resolution STEM image, (d) the magnification of a 
rectangular region with the assignment of atoms. 


800 nm thick SL (Fig. la inset) implies a smooth bound¬ 
ary between layers. Moreover, the STEM image shown in 
Fig. la proves the ability to grow high quality structures. 
In this particular case, a SL with a period of about 4 nm 
is presented. Every bright layer of the stack corresponds 
to HfNiSn, whereas dark layers demonstrate the presence 
of TiNiSn. 

Additionally, selected area electron diffraction (SAED) 
of a circular region marked in Fig. la was measured and 
the result, together with the assignment of crystalline di¬ 
rections to the reflections, is presented in Fig. lb. Every 
orange circle indicates a MgO diffraction spot, whereas 
other colors refer to film reflections. Observed diffraction 
spots confirm strongly directional growth of the film ac¬ 
cording to the orientation provided by the substrate. Due 
to the slight difference between lattice constants of both 
HH materials, we also observe a splitting of the reflection 
indicated by an arrow. Unequal lattice constants lead to 
the distortion of the lattice visible at the high-resolution 
STEM image (Fig. Ic). Based on the brightness contrast, 
one can clearly recognize piles of the constituent atoms. 

While AFM and STEM yield local information, XRD 
averages over a large area and is sensitive to inner in¬ 
terfaces, since the films are thinner than the penetration 
depth. Typical XRD data is shown in Fig. 2. The [001] 
([100]) direction of the film is parallel to the [001] ([110]) 
direction of the MgO substrate as is evident by the occur¬ 
rence of (002) and (004) film peaks. Additional satellites 
around the main diffraction lines of TiNiSn and HfNiSn 
arise from the extra periodicity present in the samples 
[17]. These satellite peaks are particularly sensitive to 
the quality of the boundaries between the constituent 
layers as any continuous fluctuations would cause broad¬ 
ening of all SL peaks resulting in difficulties to resolve 
additional oscillations [17]. Consequently, pronounced 
satellites imply better sample quality. 

In Fig. 2, besides measured XRD data, the calculation 
based on coherent x-ray scattering for a perfect SL is 
presented. To give a realistic physical shape to the simu¬ 
lated spectra, peaks were convoluted with the measured 
profiles of epitaxial TiNiSn films. Convolution param- 
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FIG. 2. Measured and calculated XRD 0 — 20 profile of a SL 
with period 21.9 nm. The intensities are rescaled to create 
an offset for clarity. Shown are the regions around (002) and 
(004) reflection. 


eters resulted from the fit of two Pseudo-Voigt profiles 
to (002) and (004) peaks for single TiNiSn films, respec¬ 
tively. Scattering factors of the atoms in the planes were 
taken from literature [18] and lattice plane spacing was 
adjusted to account for epitaxial strain. Therefore, in¬ 
stead of bulk lattice parameters of HH films (axiNiSn = 
5.941 A, ^^HfNiSn = 6.083 A) [19], the out-of-plane lat¬ 
tice parameters taken for the calculation were equal to 
<^TiNiSn = 5.906 A and anfNiSn = 6.184 A. The thick¬ 
nesses of TiNiSn and HfNiSn were set to 10.6 nm and 
11.3 nm, respectively. Comparing measured and calcu¬ 
lated spectra of the SLs, it can be concluded that not only 
the position of the peaks, but also the overall peak shape 
agree well. However, relative intensities of the peaks are 
not exactly reproduced, especially for higher order satel¬ 
lites, indicating disturbances from the perfect periodicity 
in the experimental samples. 

To measure the cross-plane thermal conductivity the 
differential 3uj method was used. The Us^; in the heater 
structure was measured by a SR850 DSP lock-in ampli- 
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FIG. 3. (a) Measured 

cross-plane thermal conduc¬ 
tivity of HH SLs with vary¬ 
ing period (black dots) and 
^calculated bascd on a series 
of thermal resistances (green 
dots). Experiment and cal¬ 
culation agree for SL periods 
in region 2. (b) Electronic 

and phononic part of the 
thermal conductivity for ev¬ 
ery SL and single films, (c) 
EWHM of the (002) rocking 
curves of the films and SLs. 


her from Stanford Research Systems. For every SL period 
two sets of regular and reference samples were prepared 
in order to confirm the reproducibility. Based on both 
pairs of results the error bars of data points were esti¬ 
mated. The difference between the thermal conductivi¬ 
ties of both series of samples having the same SL period 
was compared with an uncertainty originating from the 
measurement setup. Then the error bars were assumed 
to be equal to the greater of these values. 

Besides the measurement of the thermal conductivity, 
the XRD rocking curves (RC) of the (002) diffraction 
peak or the most intense satellite were acquired for one 
thick sample of every studied SL period. RCs were fit¬ 
ted with a Pseudo-Voigt profile and their full width at 
half maximum (FWHM) for every period was summa¬ 
rized in Fig. 3c. In contrary to our previous study on 
(Zr,Hf)NiSn/TiNiSn superlattices [7], we do not see any 
correlation between the width of RCs and the thermal 
conductivity in the TiNiSn/HfNiSn system which indi¬ 
cates the improved growth of the sytem investigated here. 

The study of room-temperature thermal conductivity 
versus the SL period is presented in Fig. 3a. The hor¬ 
izontal blue bars indicate the thermal conductivity of 
the individual film materials shown for comparison. For 
TiNiSn films we obtained an average of 4.16 W/(mK). 
Due to lower roughness and improved sample crystalline 
alignment in the current sample series this value is al¬ 
most 50% higher than the one reported by us earlier 
(2.8W/(mK)) [7]. The current result is also much closer 
to the bulk value of the thermal conductivity of TiNiSn 
reported in the literature, which ranges between (4.6 — 
9.3) W/(mK) for arc melted [2, 20, 21], 4.8 W/(mK) for 
microwave prepared [20] and 7.5W/(mK) for levitation 
melted materials [22]. Conversely, we report a reduced 
thermal conductivity down to 1.77 W/(mK) for HfNiSn 
thin films compared to literature data on bulk material 
that varies between (6.7 — 10) W/(mK) [1, 21]. As the 
mismatch between the MgO substrate and HfNiSn is con¬ 


siderably larger than for TiNiSn, the former show wider 
rocking curves, i.e. worse crystallite alignment leading to 
a reduced thermal conductance. The data shown in this 
article, however, has been gained on samples with com¬ 
parable crystalline quality. 

Black data points presented in Fig. 3a indicate that 
indeed a systematic reduction of the thermal conductivity 
was achieved. The wide range of studied SL periods may 
be divided into three regions. For large periods, labeled 
as region 3, we observe /^tot very close to a level marked 
with a red line that represents the ^Cseries of ^ bilayer 
in which the effect of internal interfaces is ignored and 
the value is calculated from the arithmetic mean of the 
measured series heat resistances of the constituent layers 


/^TiNiSn-t /^HfNiSn 

where di and d 2 are the total thicknesses of TiNiSn and 
HfNiSn. One can also notice an apparent saturation of 
the thermal conductivity for high period lengths, because 
only few interfaces (15 and 7 for periods equal to 108 nm 
and 216 nm, respectively) are not significantly affecting 
Ktot- In region 2 the reduction of the thermal conductiv¬ 
ity is much more significant and /^tot decreases systemat¬ 
ically when the number of interfaces increases from 23 to 
575. In this region, an appropriate model is the series of 
thermal resistances that include both the measured ther¬ 
mal resistivities of the individual films and the interface 
thermal resistance. To simulate the thermal conductiv¬ 
ity we apply semiclassical Boltzmann transport, based 
on the ab initio calculation of the bulk phonons using 
density-functional perturbation theory [23]. To get the 
thermal conductivity one has to integrate the phonon 
spectra over the total g-space summing up all phonon 
modes j. As an approximation we take the phonons as 
isotropic so that the angular integration gives just a con- 
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Where v{qj^x) is the group velocity of the phonon mode 
of the material x and Nq the Bose distribution function. 
It is motivated [24] that r has the form of 
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( 3 ) 


We choose Ax such that the above integral reproduces 
the experimental bulk thermal conductivity. For the in¬ 
terfaces the mean free path i = vr must be modified by 
a transmission coefficient ( and the material layer thick¬ 
ness d. This is obtained by a diffusive mismatch model 
[25] , which is based on the mismatch of the bulk phonons 
of the two materials. The frequency dependent transmis¬ 
sion through the layer and the interface is then given by 


I^X 





dNo 

dT 


dQj^x 


( 4 ) 


with an interface and thickness dependent mean free path 

^ ^ (5) 

^{Qj,x) dxC{Qj,x) AxVl^Qj^x) 

Eq. (4) gives the thermal conductivity of material x with 
layer thickness dx including the interface through which 
the phonon leaves the layer and is used in Eq. (1). The 
green circles in Eig. 3a that have been obtained from such 
modeling indicate that the interpretation of the experi¬ 
mental data using thermal interface resistances works for 
SL periods down to 3 nm (region 2 in Eig. 3a). Eor 
thinner SLs the model is no longer valid because the 
long-wavelength phonons experience an effective medium 
formed by both materials (see below), which results in a 
modified phonon dispersion as compared to each individ¬ 
ual material. 

Einally, for a SL period of 3.2 nm, the lowest ther¬ 
mal conductivity of (1.11 ±0.06) W/(mK) was measured. 
Decreasing the period even more (region 1), we observe 
an increase of /^tot due to the formation of an artificial 
crystal. Along the [100] direction the crystallographic 
unit cell of a HH compound is composed of 2 x 2 = 4 
atomic layers, where 2 neighboring atomic layers form 
the chemical unit cell. Eor the sample with the lowest 
period of 1.7 nm the new crystal consists of 1.25 unit 
cells of TiNiSn and 1.5 unit cells of HfNiSn, i.e. in total 
only 11 atomic layers. Locally the chemistry of the HH 
material requires the crystal size to be 10 or 12 atomic 
layers, whereas x-rays give the average value. In such a 
system, phonons experience the material as if it was com¬ 
posed of enlarged unit cells. In this case the still reduced 
thermal conductance must be attributed to a backfold- 
ing of the phonon dispersion relation due to the supercell 
that introduces gaps in the phonon dispersion relation at 
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FIG. 4. ZT values calculated from in-plane Seebeck coef¬ 
ficient and electrical resistivity and cross-plane thermal con¬ 
ductivity for TiNiSn, HfNiSn and their SLs. 



the new Brillouin zone boundaries and leads to an overall 
reduced dispersion of the phonon branches [26]. 

The decomposition of the thermal conductivity into 
electronic and phononic contribution is summarized in 
Fig. 3b. The electronic part of the thermal conductiv¬ 
ity /^ei, calculated according to the Wiedemann-Franz 
law, stays relatively constant for all SL periods, whereas 
the phononic part follows the trend appointed by /^tot- 
Thereby, the reduction of the cross-plane thermal con¬ 
ductivity can be definitely assigned to the phonons scat¬ 
tering at the interfaces. 

The in-plane resistivity and Seebeck coefficient were 
measured simultaneously by an LSR 3 (Linseis) in a He 
atmosphere between room temperature and 480 K. The 
plot summarizing the measured data is presented in the 
supplementary material. Due to the isoelectronic substi¬ 
tution we expect a lower interface scattering for electrons 
than for phonons and assume isotropic electronic trans¬ 
port. This assumption is confirmed by the absence of a 
systematic correlation between the electronic properties 
and the interface density. The observed increase of the 
Seebeck coefficient for very thin layers might be due to 
some degree of electron confinement in the layers. 

Using the room temperature cross-plane thermal con¬ 
ductivity and the temperature dependent in-plane power 
factor, we estimated the ZT values and present them in 
Fig. 4. The figure of merit of both single films is almost 
identical along the entire temperature range and reaches 
a value of about 0.08 at 480 K, whereas all the studied 
SLs exhibit greater ZT values. In the extreme case, for 
a SL having a period of 21.6 nm, ZT reach a value of 
almost 0.35, which is 4.4 times larger than the figure of 
merit of single films. The currently achieved value falls 
well within the interval demarcated by the previously re¬ 
ported data [4, 6, 20, 27-33] at 480 K. 

We deduce that both approaches, top-down (doping 
and phase separation) as well as bottom-up (SLs) are 
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complementary ways to enhance the figure of merit. 
While the former has no control over the inclusions size, 
the latter can vary the thickness of layers with great ac¬ 
curacy. As presented above, the figure of merit of 0.35 
for the SL is comparable to the best ZT values for bulk 
half-Heusler alloys reported in the literature at 480 K. 

In conclusion, the present research was designed to 
study the effect of interfaces on the thermal conductivity 
in HH SLs. DC sputtered TiNiSn/HfNiSn SLs reveal 
a smooth surface and well defined layered structure as 
proved by AFM, XRD and STEM. The thermal conduc¬ 
tivity of thin films was measured using the differential Scj 
method. Decreasing SL period, i.e. increasing number 
of interfaces, the expected reduction of the thermal 
conductivity is observed. For a SL period of 3.2 nm the 
minimum /^tot = (ITl ±0.06)W/(mK) was achieved. 
Further decrease of the SL period leads to the rise of /^tot 
due to the creation of an artificial crystal. We calculated 


the thermal conductivity using semiclassical Boltzmann 
transport, based on the ab initio calculation of the bulk 
phonons using density-functional perturbation theory 
and a diffusive mismatch model for phonon scattering at 
the interfaces. Experimental data and calculation agree 
for SL periods larger than 3 nm. An electron confinement 
might be the reason for the improved Seebeck coefficient 
in SLs with very low periods. The figure of merit was 
estimated based on the room temperature cross-plane 
thermal conductivity and the temperature dependent 
in-plane power factor. Due to the reduced thermal 
conductivity, the superlattices with all studied periods 
reveal enhanced ZT compared to single constituent films. 
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I. SAMPLE PREPARATION 

TiNiSn/HfNiSn superlattices (SLs) were prepared by 
DC magnetron sputtering at 520° C and an Ar pressure 
of 3.1 Pa on previously annealed MgO (001) substrates. 
While the cathodes were operating continuously at 23 W 
(TiNiSn) and 20 W (HfNiSn), a heater with the attached 
sample holder was moving back and forth between both 
cathodes to create a layered structure. The movement of 
the heater was optimized in a way that compensates a 
possible thickness gradient caused by multiple approach¬ 
ing of the cathode from the same side. Based on the 
determined deposition rate the waiting time over the re¬ 
spective target was calculated so that the thickness of 
both layers in one period was the same. 


the one dimensional heat conduction model, the thermal 
conductivity of the film is expressed by Eq. (1): 


^film — 


P • (^regular ^reference) 

26 • I • (ATi-gguiaj. ATi-gference) 


( 1 ) 


where P is the power applied to the heater, d is the thick¬ 
ness of the respective layers, 26 is the width of the heater 
line, I is the distance between the voltage probes and AT 
is represented by Eq. (2). It contains the third harmonic 
of the voltage drop 6 / 3 ^, the voltage across the heater 
line at frequency cj, i.e. the resistance of the heater 
R and and its temperature derivative. 


AT = 2 




R 

~dR 

dT 


( 2 ) 


II. 3w METHOD 

The cross-plane thermal conductivity of half-Heusler 
(HH) SLs was measured with the 3uj method [1-3]. On 
top of each SL, 200 nm of Al 2 0a^ was deposited by 
RE magnetron sputtering to electrically insulate the film 
from the heater/thermo meter line, which was produced 
in an optical lithography step from a 50 nm thick Au 
layer. Such a 20 /am wide heater line was terminated with 
two electrodes that are used to apply an AG current. Two 
additional pads were connected to the heater in between. 
An alternating current with frequency uj, applied to the 
outer electrodes, passes through the metal strip causing 
periodic Joule heating and modulating its temperature 
dependent resistance with a frequency of 2uj. Conse¬ 
quently, according to Ohm’s law, the voltage drop be¬ 
tween the inner connections contains contributions pro¬ 
portional to the first and third harmonic. The usage 
of the differential variant of the 3uj method [3] requires 
preparation of two samples - regular (1000 nm thick) and 
reference (150 nm thick) for every SL period. Assuming 
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III. THE THERMOELECTIC QUANTITIES 

The in-plane resistivity and the Seebeck coefficient 
were measured simultaneously by an LSR 3 (Linseis) in 
a He atmosphere between room temperature and 480 K. 
Both of these quantities and resulting power factor are 
summarized in Fig. 1. The horizontal arrows indicate 
the data measured for TiNiSn and HfNiSn single films 
keeping the color code presented in the legend. Due to 
the isoelectronic substitution on the M site of MNiSn we 
expect a lower interface scattering for electrons than for 
phonons and assume isotropic electronic transport. This 
assumption is confirmed by the absence of a systematic 
correlation between the electronic properties and the in¬ 
terface density. 

However, an interesting finding is an improved See¬ 
beck coefficient for SLs having extremely low periods. 
Although we used isoelectronic materials, due to the dif¬ 
ferent lattice constants a residual Fermi surface mismatch 
exists that can lead to some degree of electron confine¬ 
ment and an enhanced Seebeck coefficient. 

Having a closer look on the power factor (Fig. Ic) one 
can notice that the initial assumption about the fact that 
electronic parameters S‘^a should not be affected by the 
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FIG. 1. (a) In-plane Seebeck coefficient, (b) in-plane electrical resistivity and (c) power factor for the SLs with a given period 

and single materials TiNiSn and HfNiSn. 


SL approach is valid. None of the investigated SL samples 


exhibited lower power factor than the single HfNiSn film. 
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